program main
    implicit none
    !call DIFF1D()
    call NewDiff1D
end program